OU-TAP 161 



Effects of Formation Epoch Distribution on X-Ray Luminosity and 
Temperature Functions of Galaxy Clusters 

Motohiro Enoki Fumio Takahara ^, and Yutaka Fujita ^ 

g ■ ABSTRACT 
(N 

We investigate statistical properties of galaxy clusters in the context of hierarchical 
' clustering scenario, taking account of their formation epoch distribution, motivated by 

■ the recent finding by Fujita and Takahara that X-ray clusters form a fundamental plane, 
where the mass and formation epoch are regarded as two independent parameters. Using 

. the formalism which discriminates between major merger and accretion, the epoch of 

i a cluster formation is identified with that of the last major merger. Since tiny mass 

\ accretion after the formation does not much affect the core structure of clusters, the 

'si" ' properties of X-ray emission from clusters are determined by the total mass and density 

at their formation time. Under these assumptions, we calculate X-ray luminosity and 

■ temperature functions of galaxy clusters. We find that the behavior of luminosity 
function is different from the model which does not take account of formation epoch 
distribution, while the behavior of temperature function is not much changed. In our 

2 , model, luminosity function is shifted to a higher luminosity and shows no significant 



I evolution up to z ~ 1, independent of cosmological models. The clusters are populated 

^ ' on the temperature-luminosity plane with a finite dispersion. Since the simple scaling 

model in which the gas temperature is equal to the virial temperature fails to reproduce 
^ ' the observed luminosity-temperature relation, we also consider a model which takes the 

■ effects of preheating into account. The preheating model reproduces the observations 

much better. 

Subject headings: cosmology: theory — clusters: galaxies: general — X-rays: galaxies 



1. Introduction 

Galaxy clusters are the largest virialized objects in the universe and should be useful cosmo- 
logical probe since the properties of clusters are considered to depend on cosmological parameters 
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and to reflect the structure formation history. Galaxy clusters contain a diffuse hot gas which emits 
X-rays via thermal bremsstrahlung. The X-ray temperature and luminosity functions of clusters 
and their evolution have been used to determine the density of the universe and the amplitude of 
the rms density fluctuation. If there exist unique mass-temperature (T) and mass-luminosity (L) 
relations, one can predict X-ray temperature and luminosity functions of clusters at redshift z from 
a theoretical mass function. Therefore, it is important to investigate into correlations among phys- 
ical quantities of clusters and examine if such unique relations exist. Recently, Fujita & Takahara 
(1999a) have found that clusters at low redshifts (z < 0.1) form a plane (the fundamental plane) in 
the three dimensional space (log/9gas,0i logr^, logT), where /9gas,o is the central gas density, and 
is the core radius of clusters. The data on the plane still have a correlation and form a band on the 
plane. The observed relation L oaT^ turns out to be the cross section of the band normal to the 
major axis. The existence of the fundamental plane implies that the clusters form a two-parameter 
family, suggesting that the physical quantities of clusters are determined by the halo mass and 
density at the time of their formation and that no unique mass-temperature or mass-luminosity 
relation exists (Fujita & Takahara 1999b). In fact, using A^-body simulations, Navaro, Frcnk, & 
White (1997) found that the structure of clusters is related to their formation epoch, although they 
claimed that the clusters are a one-parameter family by assuming a unique relation between the 
mass and formation redshift. In this paper, motivated by the finding of the fundamental plane, 
we present a formulation of statistics of clusters in the context of hierarchical structure formation 
theory taking account of formation epoch distribution. 

So far. Press & Schechter (1974, hereafter PS) mass function has been widely used to calculate 
temperature and luminosity functions (e.g. Eke, Cole, & Frenk 1996). Although PS mass function 
provides the number of dark halos at a given time, it docs not contain information about the 
formation time of halos. Thus, a conventional way to compute the temperature and luminosity 
functions is simply to assume that the formation time is the same as the observed time. Since in 
this case clusters form a one-parameter family about the mass, an extension of PS theory is necessary 
to investigate the two-parameter nature of clusters, taking account of the effect of the formation 
epoch distribution. Using the merger probabilities in an extended Press & Schechter clustering 
model (Bond et al. 1991; Bower 1991), Lacey &l Cole (1993, hereafter LC) derived a formation 
epoch distribution function in an enlightening way, although LC did not calculate predictions of 
X-ray luminosity and temperature functions. The formation time in their model is defined by that 
time when the halo mass becomes half that at the observed epoch and this definition does not 
discriminate between tiny and notable relative mass capture, or between accretion and merger. 
In the hierarchical clustering scenario, low-mass objects successively merge with one another to 
build up ever more massive objects. However, major deviations from equilibrium and subsequent 
relaxation take place only when halos of comparable masses merge, while tiny mass captures have 
little effect on the capturing halos. So, it is desirable to devise such formulations. Kitayama & Suto 
(1996a,b) have attempted to describe the formation and destruction of halos within the extended PS 
prescription by discriminating between accretion and merger. They define the halo formation and 
destruction rates and define the formation epoch distribution by utilizing the survival probability. 
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Although they predict various statistical properties of X-ray clusters, they assume an empirical 
relation between temperature and luminosity in such predictions, which is not satisfactory from 
our point of view. 

In this paper, wc adopt another model of the formation epoch distribution proposed by 
Salvador-Sole, Solanes & Manrique (1999, hereafter SSM). They developed a modification of the 
extended PS model that differentiates merger from accretion and define the formation epoch as 
the epoch of the last major merger by utilizing the empirical mass-density correlation obtained by 
iV-body simulations. Following this formalism, we can obtain not only formation time distribution 
but also the halo mass at the formation time. In order to obtain statistics of galaxy clusters, we 
assume that the halo mass and density at its last major merger determine the structure of clusters. 
The temperature and luminosity of galaxy clusters are calculated in terms of the halo mass and 
density at the formation without resorting to empirical relations. On the basis of these prescrip- 
tions, we construct a simple scaling model of X-ray clusters and calculate L — T distribution, the 
temperature and luminosity functions. Since it is well known that a simple scaling model between 
gas and dark matter results in an L — T relation different from observations (Eke, Navarro, &; 
Frcnk 1998), we also investigate a simple preheating model. The paper is organized as follows. The 
formulation of statistics is presented in § 2. The X-ray cluster model is constructed in § 2.2. The 
results and discussion are described in § 3. 



2. Formulation 



2.1. Formation time distribution 



In this subsection, we summarize the SSM formalism. To follow the formation and evolution of 
halos, SSM used a modified version of the extended PS clustering model (Bond ct al. 1991; Bower 
1991, ; LC) and made a schematic distinction between minor and major mergers by defining the 
formation of a halo as the last major merger it experienced. This definition does not affect the 
abundance of halos at a given time, although it affects the description of their growth. Thus, the 
mass function is equal to the PS mass function 



n{M,t)dM = J^-P' ^^^^ 



TrMa'^{M) 



da{M) 



dM 



exp 



1 m 



dM, 



(1) 



where po is the present mean density of the universe, (5c (^) is the critical density contrast for collapse 
at t, and a{M) is the rms density fluctuation in spheres containing a mean mass M. In this paper, 
we use an approximate formula of Sc{t) for spatially fiat cosmological model (Nakamura Sz Suto 
1997) and a fitting formula of a{M) for the CDM fluctuation spectrum (Kitayama 1997). 



In the LC model, the instantaneous merger rate for halos with mass M at t per infinitesimal 
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range of final mass M' > M, or specific merger rate, is 



X exp 



dt 
(72 (M) 



1 -3/2 



da{M') 



dM' 



2 [ct2(M') C72(M)_ 



dM'. 



(2) 



In the SSM formalism, it is assumed that a halo with mass M experiences a major merger 
and is destroyed when the relative mass increment AM/M = {M' — M)/M exceeds a certain 
threshold A^. The major merger is regarded as the formation of a new halo. On the other hand, 
when AM/M < A,„,,, the event is regarded as continuous accretion; the halo keeps its identity 
and its core structure. Thus, from the specific merger rate (cq. [2]), the mass accretion rate, 
-Rmass(-^) = dM/dt, of halos with mass M at time t is defined as 



Rr, 



,(M,t)= / 

JM 



M{l+Am) 



AMr^ciM M',t)dM'. 



The destruction rate is defined as 



\M,t)= / 



M{l+Am) 



rZc{M ^ M',t)dM'. 



The formation rate should be written by 



rf[M{t),t] = +r^[M{t),t]+dMRr.^AM{t),t] 



(3) 



(4) 



(5) 



from the conservation equation for the number density of halos per unit mass along mean mass 
accretion tracks, M(i), which is the solution of the differential equation 



dM 
~dt 



i?mass[M(t),t]. 



(6) 



From this formation rate, one can obtain the distribution of formation time, tf, for halos with 

masses M at t; 



^f{tf]M,t)dtf = r^[M(t/),t/]exp|- rf[M{t'),t']drj dtf. 
Prom equation (6), the halo mass at formation time, Mj = M{tf ), becomes 

Mf = M- f R^^,[M{t'),t']dt'. 
Jtf 



(7) 



(8) 
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The value of is fixed by the fits to the empirical mass-density (or mass-radius) correlation 
obtained by A^-body simulations (Navaro, Frenk, & White 1996, 1997). The best fit is A^, = 0.6 
in a number of different cosmological models (SSM). To be precise, Mj for fixed M and Zf should 
have scatter around the mean accretion track (eq. [6]). Therefore, Mf in equation (8) is an 
approximation to the mean formation mass of true distribution. Although we ignore this scatter 
for the sake of simplicity, it could be estimated by using the algorithm given in Nusser &; Sheth 
(1999). 

In this paper, we investigate two cosmological models (SCDM and LCDM). In Table 1, we 
tabulate the cosmological density parameter (Oq), the cosmological constant (Ao), the Hubble 
constant in units of 100 km Mpc~^ (/?,), and the present rms density fluctuation in spheres of 
radius 8/i~^Mpc (cis)- Figure 1 shows the formation time distributions (eq. [7]) for the SCDM and 
LCDM models for several values of the present halo mass. Since in the LCDM model the universe 
is in the accelerated expanding stage at ^ ~ 0, the growth of density fluctuation has stopped. Thus, 
the slope of the formation time distribution in the LCDM model is less steep in comparison with 
the SCDM model at z ^ (Figure 1) . Figure 2 shows halo masses at its formation time (eq. [8]) 
for the halo with the present mass of 10^^ M©. 

Note that, in the PS formalism, there is no distinction between merger and accretion. Any 
mass capture is regarded as merger, that is, tf =t. Consequently, the formation time distribution 
is 

^f{tf,M,t) = d''itf-t), (9) 
where d^{x) is Dirac's delta function. 

2.2. X-ray cluster model 

In the SSM formalism, the mass function of dark halos at a fixed epoch is given by the PS 
mass function. On the other hand, for a given mass, the dark halos take a range of formation 
time as depicted by equation (7). Thus, dark halos form a two-parameter family. SSM found that 
the characteristic density and the scale radius of a halo at the present epoch are proportional to 
the critical density of the universe and the virial radius of the halo at the formation time t = tf, 
respectively. This means that, between major mergers, halos gradually grow through the accretion 
of surrounding matter while keeping the central part unchanged. Considering the two-parameter 
family nature of clusters, we thus assume that the physical quantities of clusters of galaxies are 
represented by the mass and the virial density at the their formation time or the last major merger. 
On this assumption, the temperature and luminosity of a cluster are determined by the formation 
redshift Zf and the halo mass Mf at zf, T = T[zf,Mf{M,Zf)], L = L[zf,Mf{M,Zf)]. Here, 
the relation between tf and Zf are straightforwardly determined once the cosmological model is 
specified. In order to obtain L and T, we use the spherical collapse model (Tomita 1969; Gunn 
Gott 1972) and assume that the cluster is spherically symmetric and intracluster gas is in an 
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isothermal hydrostatic equilibrium in the gravitational potential of dark matter halo. 



In the spherical collapse model, the virial density of a halo at the formation redshift is 
given by 



Pf = Pc{zf)Ac{zf) = pcoAcizf)- 



(10) 



where pdzf) is the critical density of the universe, Ac(z/) is the ratio of the virial density to the 
critical density, and ^{zf) is the cosmological density parameter. The index refers to the values 
at z = 0. We use the fitting formula of Bryan & Norman (1998) for the virial density of spatially 
fiat cosmological models; 



A. 



IStt^ + 82a; - 39x'^ , 



(11) 



where x = n{zf) - 1 = (J^o - l)/[nQ{l + zff + 1 - Qq]- 
The virial radius is obtained by 



(12) 



The virial temperature Tyir is given by 

IJ,mpGMf 



= 0.92 



3r/ 



2/3 



On 



n{zf)\ 



1/3 



{1+Zf) keV 



(13) 



where p is the mean molecular weight which we take to be /x = 0.6, nip is the proton mass, and kB 
is the Boltzmann constant. 



The luminosity of a cluster is 

L = 6.50 X 10-2^ 



-T— — / nAirr dr erg 
VlkeVy 7o 



(14) 



where ng is the electron number density in CGS units. Note that the gas temperature T is not 
equal to the virial temperature Tyir in general. In order to calculate luminosity from equation (14), 
we need to specify the gas density profile. Here, we adopt isothermal (3 model 



Pgas(r) = Pgas,0 



r 

1+ - 



-2.0/2 



(15) 



where Pgas,o is the central gas density and Tc is core radius. The central gas density is calculated 
through the relation 

/ pgas(r)47rr2c?r = M^^ (16) 
Jo 
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where Mgas is the total gas mass. Prom now on, we assume that Tc = r//8 (Pujita &; Takahara 
1999b). 

Thus, from (14) and (15), the luminosity is given by 

{1 + Zffflevgs-' (17) 

where fm = Mg^s/Mf is the gas mass fraction and we assume here /3 = 2/3 for simplicity. Equations 
(13) and (17) show that the temperature and luminosity are functions of the halo mass and the 
formation redshift indeed when one specifies the cosmological model and the relation between T 
and Tvir. 



2.3. Temperature and luminosity function 

Next, wc formulate statistics of galaxy clusters. Now that the temperature and luminosity are 
expressed as functions of tj and Alf, such that T = T[tf,Mf{AI,tf)] and L = L[tf,Mf{M,tf)], 
we can derive temperature and luminosity functions using the transformation of the two variables 
from {tf,M) to {T,L) as follows. 

For a mass range M ~ M + dM at t, the comoving number density of clusters formed at 
tf tf + dtf is given by 

n{tf,M;t)dtfdM = n{M;t)<^f{tf; M,t)dtfdM. (18) 

Thus the comoving number density of clusters at t with T T + dT and L L + dL is given by 

n{T,L;t)dTdL = n{tf,M;t)dtfdM. (19) 

Accordingly, 

n{T,L;t) = n{tf,M;t) §|r^ • (20) 

Since the comoving number density n{T, L; t) gives the distribution function on L — T plane, it 
reflects L — T relation of clusters. 

The temperature function and luminosity function are respectively given by 

n{T;t) = Jn{T,L;t)dL (21) 

and 

n{L;t) = J n{T,L;t)dT. (22) 

It is to be noted that these functions take account of the distribution of cluster formation time, tf. 

Thus, given the relations T = T[tf,Mf{M,tf)] and L = L[tf,Mf{M,tf)], one can calculate 
the L — T distribution function, the temperature function and the luminosity function. 



-8- 



3. Results and discussion 

In this section, we derive the temperature function, the luminosity function and L — T distri- 
bution function. We assume that the gas mass fraction of clusters fm is the same as the cosmic 
baryon ratio ^b/Clo where Jl^ is the density parameter of baryon and we adopt = 0.0125/i~^ to 
be consistent with the primordial nucleosynthesis. 

3.1. Scaling model 

First, wc examine the case where the gas temperature T is equal to the virial temperature Tyir. 
In this case, the structure is self-similar and wc call this model the scaling model. Wc calculate 
the temperature and luminosity functions based on the SSM formalism (eq. [7]) and compare 
them with those based on the PS formalism (eq. [9]). The temperature function and their redshift 
evolution are shown in Figure 3, where thin lines and symbols denote SSM and PS predictions, 
respectively. As is seen, there is little difference between the predictions of the two formalisms. 
This can be explained as follows. When wc take account of the distribution of the formation time, 
physical quantities of clusters are affected by two factors. First, if the object formed earlier, its virial 
density becomes larger (eq. [10]). Second, if the object formed earlier, its mass at the formation 
becomes smaller for a given mass at z = (Figure 2). Gas temperature depends on the virial 
density and mass at the formation time as T oc Mf/vf oc pj^'^Mj/'^, and these two factors tend 
to cancel out. This behavior of temperature function is the same as that of previous studies (e.g. 
Kitayama & Suto 1996b). This feature is common to both the cosmological models we investigate. 
The Zf dependence of temperature, T = T{zf, M = 10^^ Mq), is shown in Figure 4, which explicitly 
shows that T is roughly constant for Zf < 2. Note that Mathiesen (1999) also finds no evidence for 
a correlation between X-ray temperature and formation time in his simulation of X-ray clusters. 
The difference between SCDM model and LCDM model is mainly caused by the dependence on 
the time evolution of virial density (eq. [10]). 

The luminosity function and their redshift evolution are shown in Figure 5, where thin lines 
and symbols denote SSM and PS predictions, respectively. As is seen, SSM formalism predicts 
larger number density than PS formalism and for the SSM formalism, there is little evolution 
even in SCDM model in contrast with the PS formalism. This result is different from that of 
the temperature function, because of the difference between the dependence of temperature and 
luminosity on the virial density and mass at the formation. Since L oc pj^r^T^^"^ oc pfMfT^^"^ oc 

pj^M^J"^^ the luminosity depends more strongly on the virial density than the temperature. The 
Zf dependence of luminosity, L = L{zf,M = IO^^Mq), is presented in Figure 6, which shows 
that L increases with Zf. Therefore, earlier formed dense clusters contribute to the increase of the 
luminosity function shown in Figure 5. Moreover, the increase of L with Zf explains that there is 
little evolution of the luminosity function from z = 1.0 to z = both in SCDM and LCDM universes. 
Since the growth of density fluctuations in SCDM model is more rapid than that in LCDM model. 
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one might think that the luminosity function in SCDM model should evolve rapidly. However, 
the present result implies that it is not the case if we consider the effect of the distribution in the 
formation redshift. Thus, recent observational evidence for little evolution of luminosity function 
(Rosati et al. 1998; De Grandi et al. 1999; Nichol et al. 1999) does not necessarily mean that SCDM 
model is disfavored against LCDM model. It is also noted that if we consider the effect of formation 
epoch distribution, smaller value of erg is needed to reproduce observations compared with the PS 
formalism. 

When we compare the predictions with observations, one should note that the amplitude of 
temperature and luminosity functions can be adjusted by changing erg so that we are concerned 
with their shape. In Figure 3, the thick line is the best power-law fit to the observed low redshift 
temperature function obtained by Henry (2000). The predicted shape of the temperature function 
is consistent with observations both for SCDM and LCDM models, although agreement is better 
for LCDM model. In Figure 5 the thick line is the best fitted Schechter function to the observed 
bolometric luminosity function within z = 0.3 obtained by Ebeling et al. (1997). The predicted 
shape is much steeper than the observations both for SCDM and LCDM models. This is another 
representation of the well-known discrepancy of the L — T relation. 

Next, we investigate the distribution on the L—T relation. For several values of n[log(T), log(L); z = 
0], we plot iso-density contours on the L — T plane in Figure 7. Prom the relations T oc p'j^M^J^ 

and L oc pj^M^J"^ ^ the luminosity behaves as L oc p^j'^T'^ . As discussed in §1, in the PS formalism 
it is assumed that the observed redshift of a cluster is equal to the formation redshift (z = Zf\ 
Thus, most of the observed clusters (z ~ 0) have nearly the same virial density Pf{z) and their 
physical quantities depend only on mass. Therefore, in the PS formalism, clusters form a one- 
parameter family, which is shown by the straight line. On the other hand, in the SSM formalism, 
the virial densities take a wide range of values because depends on the formation redshift. Thus, 
clusters form a two-parameter family. The scatter of the L — T relation shown in Figure 7 reflects 
the dispersion of the halo formation time. Scharf & Mushotzky (1997) also pointed this out. Since 
Pf tends to be distributed more widely for smaller clusters, the scatter ol L — T relation is larger 
at lower T and L. Thus, the slope oi L — T relation in the SSM formalism becomes shallower than 
L ocT^. This is, however, in conflict with the observed correlation L (xT^ (e.g. David et al. 1993). 
As long as we assume that the ratio of the gas density of the core to the virial density is constant, 
this tendency persists. To resolve this discrepancy, this ratio should vary such that less massive 
clusters have smaller baryon fraction at the cluster core from which much of the X-ray emission 
originates (Metzler and Evrard 1994; Kay and Bower 1999; Valageas and Silk 1999; Wu, Fabian, 
Sz Nulsen 2000). Many authors have attributed such behavior to preheating of intracluster gas 
(Kaiser 1991; Evrard &; Henry 1991; Cavaliere et al. 1998; Balogh et al. 1999), in which intracluster 
gas had already been heated before the cluster formed. On the other hand, some authors claimed 
that intracluster gas is heated after the formation of the cluster (Loewenstein 2000; Brighenti & 
Mathews 2001). Moreover, Bryan (2000) argued that cooling of gas and the resultant galaxy 
formation are sufficient to explain to lower the gas fraction in small clusters and groups without 
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additional heating. Since these models give qualitatively similar gas distributions, we will adopt a 
preheating model in the next subsection. It is to be noted future observations may discriminate 
the heating or galaxy formation models (e.g. Fujita 2001). 



3.2. Preheating model 

We examine the effects of preheating using a simple model of Fujita & Takahara (2000) based 
on the models of Cavaliere et al. (1998). They combine effects of shock heating and preheating so 
that the gas temperature is higher than the virial temperature such that 

T = Tyi, + ^Ti, (23) 

where Ti is a given preshock temperature. In this model, (3 is given by 

3 = IliL (OA) 

^ T.ir + (3/2)ri- ^ ^ 

If Ti is not negligible compared to Tyir, /? becomes smaller, which means that hot gas expands and 
the gas density in the core decreases under the condition that the total gas mass of the cluster is 
constant. As a result, X-ray luminosity becomes smaller, which results in a steeper L — T relation. 
Comparison with observations of clusters suggests that the preheated temperature, Ti, is about 
0.5 - 2 keV (Fujita & Takahara 2000). Here, we adopt Ti = 1 keV. 



Using equation (15), we define the normalized central gas density as 

pf 

where 



rrf/rc 2 

^l(^)-4Vo (l+X^)3/^/2 ^"- (26) 

From equations (14), (15) and (25), the luminosity is written, using fd{P), as 

1/2 



L-32xl0^or V ( ) 

where we assume rf/rc = 8, and 



',2 a ^ 



{l + zfffdi/3f l2i/3)eigs-\ (27) 



^2(/3) = 4vr ^^^^^dx. (28) 

Because /? is a function of Tyir once Ti is fixed, T L are function of Zf and Mj in this 
preheating model, as is the case of the scaling model. Thus, we can calculate the Zf dependences 
of the temperature and luminosity functions using the formulation constructed in § 2.3. The Zf 
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dependence of temperature, T = T{zf,M = IO^^Mq), and that of the luminosity, L = L{zf,M = 
10"^^ Mq), are shown by thick Unes in Figures 4 and 6, respectively. Because of preheating, the 
gas temperature becomes higher and luminosity becomes lower compared with the scaling model 
predictions for this case. The temperature and luminosity functions are shown in Figures 8 and 9, 
respectively. For several values of n[log(T), log(L); z = 0], iso-density contours on the L — T plane 
are plotted in Figure 10. Because of the preheating, the gas temperature is raised and the gas 
expands compared to the case withoTit preheating. Thus, the central gas density, fd, is decreased 
and luminosity is lowered. When T^h- ^ Ti, the effect of preheating is large and luminosity is greatly 
decreased. On the other hand, when T^ir ^ Ti, a cluster is not much affected by the preheating 
and /3 ~ 1. Since we assumed /3 = 2/3 in the scaling model, the luminosity in the preheating 
model is larger than that in the scaling model in spite of preheating. Therefore, the slope of L — T 
relation in the preheating model is steeper than L oc (Figure 10), the number of clusters with 
L > lO^^crg s^^ increases (Figure 9), and the number of clusters with L < lO^^erg decreases 
(Figure 9). In the preheating model, there is also little evolution of luminosity function. This 
reason is the same as the scaling model. 

In Figures 8 and 9, the thick lines are the observed temperature and luminosity functions 
which arc the same as Figure 3 and Figure 5. Both in the SCDM and LCDM model, the slope 
of the luminosity function based on the SSM formalism are near to the observed one, although 
it is still slightly steeper in the LCDM model. On the other hand, the slope of the temperature 
function is steeper than the observed one. Better data and better preheating models are needed to 
resolve this mismatch. In Figure 10, the thick solid lines show the PS prediction which becomes 
as steep as the observed slope because of the preheating effect. The thick dotted lines represent 
the observed L — T relation obtained by David et al. (1993). Note that observed L — T relation 
has a large dispersion which is comparable to the predicted width. Both in the SCDM and LCDM 
models, the predicted L — T relation can match the observed relation rather well. 



4. Conclusions 

We have investigated the effects of formation epoch distribution on the statistical properties of 
galaxy clusters in the context of hierarchical structure formation scenario. The mass and formation 

redshift of galaxy clusters constitute two independent parameters. In this way, using the formalism 
of Salvador-Sole, Solanes & Manrique (1999), with a few plausible assumptions, we have derived 
temperature and luminosity functions and the distribution on the L—T plane. First, we investigated 
a simple scaling model in which gas temperature is equal to the virial temperature. The temperature 
function is almost the same as that in the PS formalism, while the luminosity function is shifted to 
higher luminosity and shows no significant evolution independent of the cosmological model because 
earlier formed clusters have denser intracluster gas in the cluster core. The luminosity-temperature 
relation becomes a band with a broad width instead of a linear line, but its slope becomes a little 
flatter than that of the PS formalism, which is inconsistent with the observations. Second, we have 
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examined a simple preheating model in this framework. Preheating makes the gas distributions of 

poor clusters flatter than those of rich clusters, which reduces the X-ray luminosity. The resultant 
L — T relation is steeper than that in the scaling model and becomes consistent with observations. 
Although the temperature and luminosity functions are broadly consistent with observations, too, 
better observations and better preheating models are needed for quantitative comparisons. 
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Fig. 1. — Formation redshift distribution function, M,z = 0), for several values of the present 
halo mass in (a) the SCDM model and (b) the LCDM model. 

Fig. 2. — Halo mass at the formation redshift, Mf{zf; M = 10^^ Mq), in the SCDM model and the 
LCDM model. 

Fig. 3. — Temperature function at z = and z = 1 in (a) the SCDM model and (b) the LCDM 
model. Symbols indicate the predictions based on the PS formalism, while thin lines indicate the 
predictions based on the SSM formalism. Thick line represents the power-law fit to the observed 
low redshift temperature function obtained by Henry (2000). 

Fig. 4. — The temperature-formation redshift relation: T = T{zf\M = 10^^ Mq). Thin lines are 
the relation in the scaling model, while thick lines are the relations in the preheating model. 

Fig. 5. — Luminosity function at z = and 2; = 1 in (a) the SCDM and (b) the LCDM. Symbols 
indicate the predictions based on the PS formalism, while thin lines indicate the prediction based 
on the SSM formalism. Thick line represents the best-fitting Schechter function to the observed 
bolometric luminosity function for galaxy clusters within z = 0.3 obtained by Ebeling et al. (1997). 

Fig. 6. — The luminosity-formation redshift relation: L = L{zf,M = lO^^M©). Thin lines are the 
relation in the scaling model, while thick lines are the relations in the preheating model. 

Fig. 7. — Iso-density contours on L — T plane in (a) the SCDM model and (b) the LCDM model. 
The values of log{n[log(r), log(L); 2; = 0]} are separated at equal logarithmic intervals by 1.25 
and range from —4 to —12.75. Thick line indicates the predicted L — T relation based on the PS 
formalism. This corresponds to the case Zf = and follows L oc T^. Dotted line represents the 
power-law fit to the observed L — T relation obtained by obtained by David et al. (1993). 
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Fig. 8. — Temperature function in the preheating model at z = and z = 1 in (a) the SCDM 
model and (b) the LCDM model. Thick line represents the power-law fit to the observed low 
redshift temperature function obtained by Henry (2000). 

Fig. 9. — Luminosity function in the preheating model at z = and z = 1 in (a) the SCDM model 
and (b) the LCDM model. Thick line represents the best-fitting Schechter function to the observed 
bolometric luminosity function within z = 0.3 obtained by Ebeling et al. (1997). 

Fig. 10. — Iso-density contours on the L — T plane in (a) the SCDM model and (b) the LCDM 
model. The values of log{ra[log(r), log(L); z = 0]} arc separated at equal logarithmic intervals by 
1.25 and range from —4 to —12.75. Thick line indicates the predicted L — T relation based on the 
PS formalism. This corresponds to the case of Zf = 0. Dotted line represents the power-law fit to 
the observed L — T relation obtained by obtained by David et al. (1993). 
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Table 1: Cosmological parameters of the models 
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